A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genome-wide forward genetic screening to identify receptors and proteins mediating nanoparticle uptake and intracellular processing.

Understanding how cells process nanoparticles is crucial to optimize nanomedicine efficacy. However, characterizing cellular pathways is challenging, especially if non-canonical mechanisms are involved. In this Article a genome-wide forward genetic screening based on insertional mutagenesis is applied to discover receptors and proteins involved in the intracellular accumulation (uptake and intracellular processing) of silica nanoparticles. The nanoparticles are covered by a human serum corona known to target the low-density lipoprotein receptor (LDLR). By sorting cells with reduced nanoparticle accumulation and deep sequencing after each sorting, 80 enriched genes are identified. We find that, as well as LDLR, the scavenger receptor SCARB1 also mediates nanoparticle accumulation. Additionally, heparan sulfate acts as a specific nanoparticle receptor, and its role varies depending on cell and nanoparticle type. Furthermore, some of the identified targets affect nanoparticle trafficking to the lysosomes. These results show the potential of genetic screening to characterize nanoparticle pathways. Additionally, they indicate that corona-coated nanoparticles are internalized via multiple receptors.

Alternative Splicing Is a Major Factor Shaping Transcriptome Diversity in Mild and Severe Chronic Obstructive Pulmonary Disease.

The role of alternative splicing in chronic obstructive pulmonary disease (COPD) is still largely unknown. We aimed to investigate the differences in alternatively splicing events between patients with mild-to-moderate and severe COPD compared with non-COPD control subjects and to identify splicing factors associated with aberrant alternative splicing in COPD. For this purpose, we performed genome-wide RNA-sequencing analysis of bronchial brushings from 23 patients with mild-to-moderate COPD, 121 with severe COPD, and 23 non-COPD control subjects. We found a significant difference in the frequency of alternative splicing events in patients with mild-to-moderate and severe COPD compared with non-COPD control subjects. There were from two to eight times (depending on event type) more differential alternative splicing events in the severe than in the mild-to-moderate stage. The severe COPD samples showed less intron retention and more exon skipping. It is interesting that the transcript levels of the top 10 differentially expressed splicing factors were significantly correlated with the percentage of many alternatively spliced transcripts in severe COPD. The aberrant alternative splicing in severe COPD was predicted to increase the overall protein-coding capacity of gene products. In conclusion, we observed large and significant differences in alternative splicing between bronchial samples of patients with COPD and control subjects, with more events observed in severe than in mild-to-moderate COPD. The changes in the expression of several splicing factors correlated with prevalence of alternative splicing in severe COPD. Alternative splicing can indirectly impact gene expression by changing the relative abundance of protein-coding isoforms potentially influencing pathophysiological changes. The results provide a better understanding of COPD-related alternative splicing changes.

A bronchial gene signature specific for severe COPD that is retained in the nose.

Introduction: A subset of COPD patients develops advanced disease with severe airflow obstruction, hyperinflation and extensive emphysema. We propose that the pathogenesis in these patients differs from mild-moderate COPD and is reflected by bronchial gene expression. The aim of the present study was to identify a unique bronchial epithelial gene signature for severe COPD patients. Methods: We obtained RNA sequencing data from bronchial brushes from 123 ex-smokers with severe COPD, 23 with mild-moderate COPD and 23 non-COPD controls. We identified genes specific to severe COPD by comparing severe COPD to non-COPD controls, followed by removing genes that were also differentially expressed between mild-moderate COPD and non-COPD controls. Next, we performed a pathway analysis on these genes and evaluated whether this signature is retained in matched nasal brushings. Results: We identified 219 genes uniquely differentially expressed in severe COPD. Interaction network analysis identified VEGFA and FN1 as the key genes with the most interactions. Genes were involved in extracellular matrix regulation, collagen binding and the immune response. Of interest were 10 genes (VEGFA, DCN, SPARC, COL6A2, MGP, CYR61, ANXA6, LGALS1, C1QA and C1QB) directly connected to fibronectin 1 (FN1). Most of these genes were lower expressed in severe COPD and showed the same effect in nasal brushings. Conclusions: We found a unique severe COPD bronchial gene signature with key roles for VEGFA and FN1, which was retained in the upper airways. This supports the hypothesis that severe COPD, at least partly, comprises a different pathology and supports the potential for biomarker development based on nasal brushes in COPD.

FastCAR: fast correction for ambient RNA to facilitate differential gene expression analysis in single-cell RNA-sequencing datasets.

Cell type-specific differential gene expression analyses based on single-cell transcriptome datasets are sensitive to the presence of cell-free mRNA in the droplets containing single cells. This so-called ambient RNA contamination may differ between samples obtained from patients and healthy controls. Current ambient RNA correction methods were not developed specifically for single-cell differential gene expression (sc-DGE) analyses and might therefore not sufficiently correct for ambient RNA-derived signals. Here, we show that ambient RNA levels are highly sample-specific. We found that without ambient RNA correction, sc-DGE analyses erroneously identify transcripts originating from ambient RNA as cell type-specific disease-associated genes. We therefore developed a computationally lean and intuitive correction method, Fast Correction for Ambient RNA (FastCAR), optimized for sc-DGE analysis of scRNA-Seq datasets generated by droplet-based methods including the 10XGenomics Chromium platform. FastCAR uses the profile of transcripts observed in libraries that likely represent empty droplets to determine the level of ambient RNA in each individual sample, and then corrects for these ambient RNA gene expression values. FastCAR can be applied as part of the data pre-processing and QC in sc-DGE workflows comparing scRNA-Seq data in a health versus disease experimental design. We compared FastCAR with two methods previously developed to remove ambient RNA, SoupX and CellBender. All three methods identified additional genes in sc-DGE analyses that were not identified in the absence of ambient RNA correction. However, we show that FastCAR performs better at correcting gene expression values attributed to ambient RNA, resulting in a lower frequency of false-positive observations. Moreover, the use of FastCAR in a sc-DGE workflow increases the cell-type specificity of sc-DGE analyses across disease conditions.

FIRRM/C1orf112 is synthetic lethal with PICH and mediates RAD51 dynamics.

Joint DNA molecules are natural byproducts of DNA replication and repair. Persistent joint molecules give rise to ultrafine DNA bridges (UFBs) in mitosis, compromising sister chromatid separation. The DNA translocase PICH (ERCC6L) has a central role in UFB resolution. A genome-wide loss-of-function screen is performed to identify the genetic context of PICH dependency. In addition to genes involved in DNA condensation, centromere stability, and DNA-damage repair, we identify FIGNL1-interacting regulator of recombination and mitosis (FIRRM), formerly known as C1orf112. We find that FIRRM interacts with and stabilizes the AAA+ ATPase FIGNL1. Inactivation of either FIRRM or FIGNL1 results in UFB formation, prolonged accumulation of RAD51 at nuclear foci, and impaired replication fork dynamics and consequently impairs genome maintenance. Combined, our data suggest that inactivation of FIRRM and FIGNL1 dysregulates RAD51 dynamics at replication forks, resulting in persistent DNA lesions and a dependency on PICH to preserve cell viability.

Transcriptomic Comparison of Human Peripartum and Dilated Cardiomyopathy Identifies Differences in Key Disease Pathways.

Peripartum cardiomyopathy (PPCM) is a rare form of acute onset heart failure that presents in otherwise healthy pregnant women around the time of delivery. While most of these women respond to early intervention, about 20% progress to end-stage heart failure that symptomatically resembles dilated cardiomyopathy (DCM). In this study, we examined two independent RNAseq datasets from the left ventricle of end-stage PPCM patients and compared gene expression profiles to female DCM and non-failing donors. Differential gene expression, enrichment analysis and cellular deconvolution were performed to identify key processes in disease pathology. PPCM and DCM display similar enrichment in metabolic pathways and extracellular matrix remodeling suggesting these are similar processes across end-stage systolic heart failure. Genes involved in golgi vesicles biogenesis and budding were enriched in PPCM left ventricles compared to healthy donors but were not found in DCM. Furthermore, changes in immune cell populations are evident in PPCM but to a lesser extent compared to DCM, where the latter is associated with pronounced pro-inflammatory and cytotoxic T cell activity. This study reveals several pathways that are common to end-stage heart failure but also identifies potential targets of disease that may be unique to PPCM and DCM.

Smoking increases expression of the SARS-CoV-2 spike protein-binding long ACE2 isoform in bronchial epithelium.

After more than two years the COVID-19 pandemic, that is caused by infection with the respiratory SARS-CoV-2 virus, is still ongoing. The risk to develop severe COVID-19 upon SARS-CoV-2 infection is increased in individuals with a high age, high body mass index, and who are smoking. The SARS-CoV-2 virus infects cells of the upper respiratory tract by entering these cells upon binding to the Angiotensin-converting enzyme 2 (ACE2) receptor. ACE2 is expressed in various cell types in the lung but the expression is especially high in goblet and ciliated cells. Recently, it was shown that next to its full-length isoform, ACE2 also has a short isoform. The short isoform is unable to bind SARS-CoV-2 and does not facilitate viral entry. In the current study we investigated whether active cigarette smoking increases the expression of the long or the short ACE2 isoform. We showed that in active smokers the expression of the long, active isoform, but not the short isoform of ACE2 is higher compared to never smokers. Additionally, it was shown that the expression of especially the long, active isoform of ACE2 was associated with secretory, club and goblet epithelial cells. This study increases our understanding of why current smokers are more susceptible to SARS-CoV-2 infection, in addition to the already established increased risk to develop severe COVID-19.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Liver transcriptomic and methylomic analyses identify transcriptional mitogen-activated protein kinase regulation in facultative hibernation of Syrian hamster.

Hibernation consists of alternating torpor-arousal phases, during which animals cope with repetitive hypothermia and ischaemia-reperfusion. Due to limited transcriptomic and methylomic information for facultative hibernators, we here conducted RNA and whole-genome bisulfide sequencing in liver of hibernating Syrian hamster (Mesocricetus auratus). Gene ontology analysis was performed on 844 differentially expressed genes and confirmed the shift in metabolic fuel utilization, inhibition of RNA transcription and cell cycle regulation as found in seasonal hibernators. Additionally, we showed a so far unreported suppression of mitogen-activated protein kinase (MAPK) and protein phosphatase 1 pathways during torpor. Notably, hibernating hamsters showed upregulation of MAPK inhibitors (dual-specificity phosphatases and sproutys) and reduced levels of MAPK-induced transcription factors (TFs). Promoter methylation was found to modulate the expression of genes targeted by these TFs. In conclusion, we document gene regulation between hibernation phases, which may aid the identification of pathways and targets to prevent organ damage in transplantation or ischaemia-reperfusion.

Airway wall splice quantitative trait locus analysis reveals novel downstream mechanisms for known asthma single-nucleotide polymorphisms.

Studying the effects of asthma SNPs on alternative splicing can lead to new insights into asthma pathophysiology. More specifically, a 17q12 SNP is associated to alternative splicing of GSDMB. https://bit.ly/3W49oTs.

A proteomics approach to identify COPD-related changes in lung fibroblasts.

Lung fibroblasts are implicated in abnormal tissue repair in chronic obstructive pulmonary disease (COPD). The exact mechanisms are unknown and comprehensive analysis comparing COPD- and control fibroblasts is lacking. The aim of this study is to gain insight into the role of lung fibroblasts in COPD pathology using unbiased proteomic and transcriptomic analysis. Protein and RNA were isolated from cultured parenchymal lung fibroblasts of 17 patients with stage IV COPD and 16 non-COPD controls. Proteins were analyzed using LC-MS/MS and RNA through RNA sequencing. Differential protein and gene expression in COPD was assessed via linear regression, followed by pathway enrichment, correlation analysis, and immunohistological staining in lung tissue. Proteomic and transcriptomic data were compared to investigate the overlap and correlation between both levels of data. We identified 40 differentially expressed (DE) proteins and zero DE genes between COPD and control fibroblasts. The most significant DE proteins were HNRNPA2B1 and FHL1. Thirteen of the 40 proteins were previously associated with COPD, including FHL1 and GSTP1. Six of the 40 proteins were related to telomere maintenance pathways, and were positively correlated with the senescence marker LMNB1. No significant correlation between gene and protein expression was observed for the 40 proteins. We hereby describe 40 DE proteins in COPD fibroblasts including previously described COPD proteins (FHL1, GSTP1) and new COPD research targets like HNRNPA2B1. Lack of overlap and correlation between gene and protein data supports the use of unbiased proteomics analysis and indicates that different types of information are generated with both methods.

Sister chromatid exchanges induced by perturbed replication can form independently of BRCA1, BRCA2 and RAD51.

Sister chromatid exchanges (SCEs) are products of joint DNA molecule resolution, and are considered to form through homologous recombination (HR). Indeed, SCE induction upon irradiation requires the canonical HR factors BRCA1, BRCA2 and RAD51. In contrast, replication-blocking agents, including PARP inhibitors, induce SCEs independently of BRCA1, BRCA2 and RAD51. PARP inhibitor-induced SCEs are enriched at difficult-to-replicate genomic regions, including common fragile sites (CFSs). PARP inhibitor-induced replication lesions are transmitted into mitosis, suggesting that SCEs can originate from mitotic processing of under-replicated DNA. Proteomics analysis reveals mitotic recruitment of DNA polymerase theta (POLQ) to synthetic DNA ends. POLQ inactivation results in reduced SCE numbers and severe chromosome fragmentation upon PARP inhibition in HR-deficient cells. Accordingly, analysis of CFSs in cancer genomes reveals frequent allelic deletions, flanked by signatures of POLQ-mediated repair. Combined, we show PARP inhibition generates under-replicated DNA, which is processed into SCEs during mitosis, independently of canonical HR factors.

A survey of current methods to detect and genotype inversions.

Polymorphic inversions are ubiquitous in humans and they have been linked to both adaptation and disease. Following their discovery in Drosophila more than a century ago, inversions have proved to be more elusive than other structural variants. A wide variety of methods for the detection and genotyping of inversions have recently been developed: multiple techniques based on selective amplification by PCR, short- and long-read sequencing approaches, principal component analysis of small variant haplotypes, template strand sequencing, optical mapping, and various genome assembly methods. Many methods apply complex wet lab protocols or increasingly refined bioinformatic analyses. This review is an attempt to provide a practical summary and comparison of the methods that are in current use, with a focus on metrics such as the maximum size of segmental duplications at inversion breakpoints that each method can tolerate, the size range of inversions that they recover, their throughput, and whether the locations of putative inversions must be known beforehand.

Chromosomal Instability Characterizes Pediatric Medulloblastoma but Is Not Tolerated in the Developing Cerebellum.

Medulloblastoma is a pediatric brain malignancy that consists of four transcriptional subgroups. Structural and numerical aneuploidy are common in all subgroups, although they are particularly profound in Group 3 and Group 4 medulloblastoma and in a subtype of SHH medulloblastoma termed SHHα. This suggests that chromosomal instability (CIN), the process leading to aneuploidy, is an important player in medulloblastoma pathophysiology. However, it is not known if there is ongoing CIN in medulloblastoma or if CIN affects the developing cerebellum and promotes tumor formation. To investigate this, we performed karyotyping of single medulloblastoma cells and demonstrated the presence of distinct tumor cell clones harboring unique copy number alterations, which is suggestive of ongoing CIN. We also found enrichment for processes related to DNA replication, repair, and mitosis in both SHH medulloblastoma and in the highly proliferative compartment of the presumed tumor cell lineage-of-origin, the latter also being sensitive to genotoxic stress. However, when challenging these tumor cells-of-origin with genetic lesions inducing CIN using transgenic mouse modeling, we found no evidence for large chromosomal aberrations in the cerebellum or for medulloblastoma formation. We therefore conclude that without a background of specific genetic mutations, CIN is not tolerated in the developing cerebellum in vivo and, thus, by itself is not sufficient to initiate medulloblastoma.

GeneNetTools: tests for Gaussian graphical models with shrinkage.

MOTIVATION: Gaussian graphical models (GGMs) are network representations of random variables (as nodes) and their partial correlations (as edges). GGMs overcome the challenges of high-dimensional data analysis by using shrinkage methodologies. Therefore, they have become useful to reconstruct gene regulatory networks from gene-expression profiles. However, it is often ignored that the partial correlations are 'shrunk' and that they cannot be compared/assessed directly. Therefore, accurate (differential) network analyses need to account for the number of variables, the sample size, and also the shrinkage value, otherwise, the analysis and its biological interpretation would turn biased. To date, there are no appropriate methods to account for these factors and address these issues. RESULTS: We derive the statistical properties of the partial correlation obtained with the Ledoit-Wolf shrinkage. Our result provides a toolbox for (differential) network analyses as (i) confidence intervals, (ii) a test for zero partial correlation (null-effects) and (iii) a test to compare partial correlations. Our novel (parametric) methods account for the number of variables, the sample size and the shrinkage values. Additionally, they are computationally fast, simple to implement and require only basic statistical knowledge. Our simulations show that the novel tests perform better than DiffNetFDR-a recently published alternative-in terms of the trade-off between true and false positives. The methods are demonstrated on synthetic data and two gene-expression datasets from Escherichia coli and Mus musculus. AVAILABILITY AND IMPLEMENTATION: The R package with the methods and the R script with the analysis are available in https://github.com/V-Bernal/GeneNetTools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The developmental stage of the medulloblastoma cell-of-origin restricts Sonic hedgehog pathway usage and drug sensitivity.

Sonic hedgehog (SHH) medulloblastoma originates from the cerebellar granule neuron progenitor (CGNP) lineage, which depends on Hedgehog signaling for its perinatal expansion. Whereas SHH tumors exhibit overall deregulation of this pathway, they also show patient age-specific aberrations. To investigate whether the developmental stage of the CGNP can account for these age-specific lesions, we analyzed developing murine CGNP transcriptomes and observed highly dynamic gene expression as a function of age. Cross-species comparison with human SHH medulloblastoma showed partial maintenance of these expression patterns, and highlighted low primary cilium expression as hallmark of infant medulloblastoma and early embryonic CGNPs. This coincided with reduced responsiveness to upstream SHH pathway component Smoothened, whereas sensitivity to downstream components SUFU and GLI family proteins was retained. Together, these findings can explain the preference for SUFU mutations in infant medulloblastoma and suggest that drugs targeting the downstream SHH pathway will be most appropriate for infant patients.